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Abstract 

We present calculations of the complex-temperature zeros of the partition functions for 2D 
Ising models on the square lattice with spin s = 1, 3/2, and 2. These give insight into 
complex-temperature phase diagrams of these models in the thermodynamic limit. Support 
is adduced for a conjecture that all divergences of the magnetisation occur at endpoints of 
arcs of zeros protruding into the FM phase. We conjecture that there are 4[s 2 ] — 2 such arcs 
for s > 1, where [x] denotes the integral part of x. 
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The Ising model has long served as a simple prototype of a statistical mechanical sys- 
tem which (for spatial dimensionality d > 1) undergoes a phase transition with associated 
spontaneous symmetry breaking (SSB) and long range order. The zero-field d = 2 spin 1/2 
Ising model is exactly solvable; the free energy / and spontaneous magnetisation M were 
first derived by Onsager and Yang, respectively Jj], || (both for the square lattice). However, 
no exact closed-form expressions have ever been found for / or any other thermodynamic 
quantity for any higher-spin Ising model in d — 2 (or higher d). It is therefore of con- 
siderable value to establish further properties of higher-spin Ising models, since these help 
one to make progress toward an exact solution. In this letter we present calculations of the 
complex-temperature (CT) zeros of the partition functions of the 2D Ising model on the 
square lattice for the higher spin values s = 1, 3/2, and 2. These results enable one to infer 
information about the CT phase diagrams in the thermodynamic limit for these models. 
Combining our results with analyses of low-temperature series by Guttmann et al. ||, |J , we 
adduce support for our conjecture on CT divergences in M. We also infer a new conjecture 
on arcs of zeros. 

There are several reasons for interest in CT properties of spin models: (i) one can under- 
stand more deeply the behaviour of various thermodynamic quantities by seeing how they 
behave as analytic functions of CT (e.g., CT singularities can affect behaviour for physical 
temperature); (ii) one can see how the physical phases of a given model generalise to re- 
gions for CT; (iii) as noted, a knowledge of the CT properties helps in the search for exact 
solutions. The earliest papers on CT properties in the Ising model include Refs. ||-[@ for 
s = 1/2; CT singularities for higher-spin (2D and 3D) Ising models were first studied in Ref. 
H (see also Refs. |I| PI). Later papers for s = 1/2 include Refs. Hl2|-||l7|. 



The spin s (nearest-neighbour) Ising model on the square lattice is defined, in standard 
notation, by the partition function Z = J2{s n } ^^ n , with 

H = -(J/s 2 ) ]T S n S n ,-(H/s)J2Sn (1) 

<nn'> n 

where S n G {— s, —s + 1, ...,s — l,s} and (3 = (/c^T) -1 . H = unless otherwise indicated. 
We define K = (3 J, v = tannic, and Us = e~ K ' s2 . Z is then a generalised (i.e. with 
negative as well as positive powers) polynomial in u s . We denote the physical critical point 
as K c for a given s, and (u s ) c = exp(—K c /s 2 ). The (reduced) free energy is / = — (3F = 
limjv^oo iV" 1 In Z in the thermodynamic limit. 

Early series studies (mainly for 3D) gave evidence that the critical exponents of the phase 
transition between the ^-symmetric, paramagnetic (PM) phase and the phase with SSB 
and long range ferromagnetic (FM) order are independent of s 110, IB]. This is in agreement 
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with expectations from renormalisation group arguments, since changing s does not change 
the Z 2 symmetry group of 7i. However, for CT, the value of s does have interesting effects. 

Since there is no closed-form expression for / for s > 1, one cannot directly determine 
the locus of CT points where / is non-analytic (aside from the points K = ±00 where it is 
trivially non-analytic). However, from calculations of zeros for the 2D spin 1/2 Ising model 
and comparison with exact results, one knows that as the lattice size increases, the zeros of 
Z occur on, or progressively closer to, the above locus of points where, in the thermodynamic 
limit, / is non-analytic. By studying the zeros of Z on finite lattices with varying boundary 
conditions (BC), one can thus make reasonable inferences about this locus of points and the 
corresponding CT phase diagram. 

For the d = 2, s = 1/2 Ising model with isotropic couplings J on the square lattice (as 
well as on the triangular and honeycomb lattices, and also certain heteropolygonal lattices 
|17j)> as N s — > 00, these zeros of the partition function in z merge to form continuous curves, 
including possible line segments. The CT phase diagram is comprised of phases bounded by 
such curves where leading and next-to-leading eigenvalues of the transfer matrix T (which 
are, in general, complex here) become degenerate in magnitude and hence there is a non- 
analytic change in /. We show now that (for isotropic J) the zeros of Z again merge to form 
a one- dimensional variety (i.e., curves including possible line segments, instead of areas) in 
the complex u s plane for arbitrary s. Using the the discrete translation invariance of the 
square lattice, one can, by Fourier transformation methods, write 

where A, B are functions of u s . The locus of points where g and hence / are non-analytic 
is given by an equation involving only the argument of g, namely, A + 5(cos^i + cos #2)- 
Now since the two independent (periodic) variables 6± and 82 only enter in the combination 
cos#i + cos#2 = x, the argument of g can be written as A(u s ) + B(u s )x, so that the equation 
involves a single independent real variable (—2 < x < 2). Hence, the continuous locus of 
points where / is non-analytic (i.e., excluding the trivial isolated isolated points K = ±00) 
is a one-dimensional variety in the u s plane. The zeros of Z and corresponding curves where 
/ is non-analytic are invariant under u s —>■ u* and also, for a bipartite lattice, under the 
mapping u s — > l/u s . 

In order to compute the zeros, we calculate Z for finite lattices with specified BC's. We 



have done this by means of a transfer matrix method [113 , and have used both periodic and 
helical BC's (resp. PBC, HBC). We have also used a type of helical BC introduced by Creutz 
(CHBC) which reduces finite-size effects. We follow the notation of Ref. |2D| to label 
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the lattice in this case: N\ x (N s /Ni), where N s denotes the total number of sites and Ni 
denotes the number of sites in all the rows except the last. 

In Fig. 1 we show the zeros of Z for the Ising model on a square lattice with (i) s = 1, 
lattice size 8 x 10, PBC; (ii) s = 3/2, lattice size 6x6, PBC; (iii) s = 2, lattice 5 x 
(38/5), CHBC The horizontal and vertical axes are Re(u s ) and Im(u s ) for each s. We have 
calculated zeros for a number of different lattices. Although we only show one plot for each 
s here, the features which we point out were observed on all of the lattices used, with both 
PBC's and (C)HBC's. 

First, up to some slight scatter, the zeros can be connected by curves, in agreement with 
our general argument above. Second, for each s, the model has physical PM, FM, and AFM 
phases, and these have analytic continuations to form CT phases in the u s plane; we denote 
these as (CT) PM, FM, and AFM phases. The symmetry of the phases about the unit circle 
\u s \ — 1 follows from the u s — > l/u s symmetry of the model. 

Third, the CT phase diagram also includes phases (labelled "O" ) which have no overlap 
with any physical phase. For s = 1,2 there is an Oi phase to the left of the FM phase 
separating it from the AFM phase. We conjecture that this is true for abitrary integral s. 
We restrict this conjecture to integral s since we already know that it does not hold in the 
exactly solved s = 1/2 model; in that case, there is no O phase (in the U\/2 = e~ 4K variable), 
and the FM and AFM (and PM) phases are directly contiguous at U1/2 = —1 (see Fig. 1(c) 
of Ref. Jlij]). The number of O phases increases as a function of s. For s = 1, our results 
suggest that in addition to Oi, there is also an O2 above, and its complex-conjugate (c.c.) 

phase below, the FM phase. There is no definite indication of any other O phases. For 
s = 3/2, there is an Oi phase or a c.c. pair of Oi and phases to the left of the FM phase, 
and c.c. O2 and phases to the upper and lower left of the FM phase. There may be other 
O phases also. For s = 2, there appear to be several c.c. pairs of O phases in addition to 
the definite Oi phase. 

Fourth, one may study special points separating CT phases. Because of the u s — > l/u s 
symmetry, and the consequent manner in which the phases group themselves around the 
unit circle \u s \ = 1, it follows that if two such phases are contiguous at a single point, then 
this point lies on this unit circle. An exactly known example is the case U1/2 — — 1 for 
s = 1/2. From our s = 1 zeros, it is likely that one such point u^u separates the Oi and 
O2 phase (so that also u*. 12 separates the Oi and phases); if this is the case, then it is 
likely that u\-\2 = e 27ri//3 . For s = 3/2, we infer likely intersection points of phase boundaries 
at u 3/2 -a = i and u 3/ 2; b = e 4 "/ 5 = (1/4) (-V5 - 1 + iy/2^h- y/Z) ~ -0.809 + 0.588i 
together with their c.c.'s. For s = 2, we infer likely intersection points at U2- a = e 7 ™/ 3 , 
u 2 . b = e 2 ^/ 5 = (1/4) ( -1 + V5 + zV2^5 + ~ 0.309 + 0.951z, and w 2;c = e 4wi/5 , together 
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with their c.c's. There may be other such points for s = 2. 

Fifth, in addition to the zeros that lie on boundaries which, in the thermodynamic limit, 
would separate CT phases, there are also zeros which lie on arcs which do not separate 
phases, but rather protrude into, and end in, the FM and AFM phases. As our exact 
solution of the CT phase diagram for the ID spin s Ising model shows, such arcs can 
represent the degeneracy, on a finite curve, of the leading and next-to-leading eigenvalues 
of T. For s — 1/2, there are no such protruding arcs for the square lattice, although exact 
results show that they do occur for other 2D lattices (triangular (in FM), honeycomb (PM), 
kagome (PM), 3-12 lattice (PM) [IT]])- We denote the number of arcs protruding into the 
phase ph (and terminating in endpoints (e)) as N e)P h. If these arcs occur at complex values 
of u s , then the u s — > u* symmetry implies that they must occur as c.c. pairs, so N £tFM is 
even. For all of the s values that we have studied (on the square lattice), the arcs do occur 
at complex values of u s , which leads one to infer that N e pM is even on the square for all s 
(this will follow from our conjecture Independent of this, for any bipartite lattice the 
u s — > l/u s symmetry implies that N Sj fm = N b> afm- 

For s — 1 (Fig. 1(a)), we find an arc protruding in a "NW-SE" direction into the 
FM phase and a corresponding arc protruding in a SE-NW direction into the AFM phase, 
together with their complex conjugates, so N e>FM = 2 for s = 1. For s = 3/2 (Fig. 1(b)), 
we find three c.c. pairs of arcs protruding into the FM phase, so N Bi fm = 6. For s = 2 (Fig. 
1(c)), we infer seven such c.c. pairs of arcs, so N £jFM = 14. Our results are summarised in 
Table 1.0 For s > 3/2 we list the positions of the last zeros to an accuracy of 0.01 in real 
and imaginary parts because this is typical of the differences between the locations of these 
zeros for different lattice sizes and BC's. The values of |(M s ) e ,z| are calculated using these 
positions to their full accuracy for the lattices in Fig. 1. The values of (u s ) c used for the last 
column are from Refs. |3], f|]. 

From previous studies of the s — 1/2 Ising model for lattices in d = 2,3, we have been 
led to formulate several general conjectures for CT singularities in M and the susceptibility 
X- One of these is Conj. 1: For the zero-field s — 1/2 Ising model on an arbitrary (regular) 
lattice with even coordination number q, M diverges at a point z = e~ 2K if and only if z is 
an endpoint of an arc (of non-analyticities of /, as discussed above) protruding into the FM 
phase. A corollary of Conj. 1 is that on a bipartite lattice, for each such endpoint there is 
another with K — > —K at which the staggered magnetisation M st diverges.^] In Ref. []TB 



we 



1 This inference does not apply to other lattices; indeed, one knows of cases where N ei FM is odd: e.g., 
exact results show that N e< FM — 1 for s — 1/2 on the triangular lattice. 

2 We add the usual caution that analyses on finite lattices up to a given size can never rigorously exclude 
the possibility that there are changes in the arcs as one gets closer to the thermodynamic limit. 

3 We record a related conjecture here Il7): Conj. 2: For the s = 1/2 Ising model on an arbitrary (regular) 
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s 


N e ,FM 


(Ua)e,Z 


\{Us)e,z\ 
(«s)c 


1/2 









1 


2 


-0.305 ±0.381z 


0.874 


3/2 


6 


0.09±0.65i 
-0.07±0.72i 
-0.54 ± 0.34i 


0.88 
0.98 
0.85 


2 


14 


0.38 ± 0.65i 
0.30±0.73i 
0.21 ±0.80z 
-0.23±0.69i 
-0.40 ± 0.69* 
-0.68 ± 0.49i 
-0.655 ± 0.29i 


0.90 
0.95 

1.0 
0.87 
0.965 

1.0 
0.86 



Table 1: Endpoints of arcs protruding into the FM phase. Entries in column denoted (w s ) e ,z 
are the positions of the last zeros on the arcs in Figs. l(a)-(c). 

proved (as Theorem 1 in the first paper) that if M diverges at a given point z = e~ 2K , so 
does Xi an d thus, for a bipartite lattice, if M st diverges, so does the staggered susceptibility, 
X^-Q Since there is no obvious reason why this conjecture should be limited to s — 1/2, it 
is natural to consider the following generalisation: Conj. Is: For the zero- field Ising model 
with arbitrary s on an arbitrary (regular) lattice with even q, M diverges at a point u s if and 
only if u s is an endpoint of an arc of non-analyticities of / protruding into the FM phase. 
A corollary of Conj. Is is that on a bipartite lattice, to each such point, there corresponds 
a point l/u s in the AFM phase, and M st diverges at a point l/u s in the AFM phase if and 
only if l/u s is the endpoint of an arc protruding into the AFM phase (the image, under 
u s — > l/u s , of the arc in the FM phase). 

Combining our determination of the positions of the endpoints of arcs protruding into 
the FM phase with the results of low-temperature series analyses in Refs. || (s = 1) and 
H 1 < s < 3, we find that for the values s = 1, 3/2, and 2 where the comparison can be 
made, there is excellent support for Conjecture Is. Of course, the position of the last zero 
on an arc calculated for a finite lattice will not, in general, be precisely the same as the 



lattice with odd coordination number q, in addition to the divergences at arc endpoints specified by Conj. 
1, M also diverges at the point z = — 1 if and only if this point can be reached from within the FM phase. 



4 Our study of the exactly solvable s = 1/2 Ising model with (3H = in/2 in Ref. |15 suggests that it is 
possible to generalise Conj. 1 to a class of nonzero external fields. 
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endpoint in the thermodynamic limit, but we expect it to be close for the lattices of the sizes 
that we are using. For s = 1, Ref. || found CT divergences in M at the c.c points U\ = 
—0.30196(3) ±0.37875(5)2, which match nicely with the arc endpoints we have found for this 
case. As a consequence of our Theorem 1 in Ref. these divergences in M automatically 



imply divergences in x, and these were also observed in Ref. ||. For s = 3/2 and 2, all of 
the divergences in M again match the locations of the arc endpoints which we have found 
0. For example, for s = 3/2, preliminary results give [|J u 3 / 2 = 0.0948(1) ± 0.6410(2)z and 
—0.5291(2) ±0.3379(2)i for the two such singularities closest to the origin, in agreement with 
our values in Table 1, with a similar match for the third. These comparisons indicate that 
the last zero on an arc for the finite lattices which we have used is a distance of about 0.01 
farther away from the origin than the singularity positions from series analyses 0, |3|], which 
presumably reflect the thermodynamic limit. (For s — 1, the difference in distances is less, 
about 0.004.) We infer that the last zeros on the arcs give values of p = \u s \/{u s ) c which are 
larger than the exact values by roughly 1 %. 

These results differ with an old conjecture || that for the spin s Ising model the number 
N cs ,p<i of complex-temperature singularities within the disk p < 1 is N CSjP< i = qs — 2, where 
q is the coordination number of the lattice. Here, this would yield N CSjP< i = 4s — 2 = 2, 4, 
and 6 for s = 1, 3/2, and 2, rather than the respective values 2, 6, and > 10 which we find 
with p < 1. 

Combining our present calculation of partition function zeros with an exact determination 
of the CT phase structure of the ID spin s Ising model [pi], we make the following conjecture: 
For the spin s Ising model on the square lattice, the number of endpoints of arcs protruding 
into the CT FM (equivalently, CT AFM phase) is 

N e , FM = N eAFM = 4[s 2 ] - 2 (3) 

for s > 1, where [x] denotes the integral part of x. Combining this with the exactly known 
result for s = 1/2, the general formula reads N £tFM = N eAFM = max{ 0,4[s 2 ] -2}. 

We next discuss the s dependence of the PM phase. In general, for a fixed lattice 
and temperature, an increase in s has the effect of allowing more randomness, and hence 
increasing the disorder in the model. This is reflected in the well-known fact that K c is a 
monotonically increasing [e~ Kc is a monotonically decreasing) function of s ||, 10, 18), |J. 



Thus, the disordered, PM phase increases in extent as s increases. However, although K c 
increases as a function of s, the ratio K c /s 2 which occurs in the variable u s decreases, 
approaching zero as s — > oo. Indeed, as s -> oo, K c approaches a finite limiting value, 
K c (s) y K c (s = oo). This follows since in this s — > oo limit (renormalising Z so that each 
spin summation is normalised to measure 1 rather than 2s + 1), one obtains the continuous- 
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spin Ising model, with Z = (J[ n J_ 1 (ds n /2)) exp(K J2<nn'> ^nS n ')- For d > 1 this model has 
a transition at a finite value of K |[22|| , viz., lim^oo K c (s). Since K c /s 2 — > const. / s 2 — > as 
s — > oo, it follows that (ii s ) c / 1 as s -> oo. Hence, in terms of tt s , the physical PM phase 
appears to decrease in size as s increases and, in the limit as s —>■ oo, to contract to a set of 
measure zero. This has also been noted in Ref. |§], where it was suggested that as s — > oo, 
both the physical critical point and the CT singularities in u s will converge to the unit circle 
\u s \ = 1. In the case of the physical critical point, one should note, however, that this is 
really due to the definition of u s ; in fact, in a fixed variable such as K or v, the physical PM 
phase actually increases in extent, reaching a limiting size, as s — > oo. From our plots of 
zeros, we conclude more generally that not just the PM phase but the O phases also cluster 
more closely around the unit circle \u s \ = 1 as s increases. We can thus state a stronger 
form of the conjecture of Ref. ||: As s —>■ oo, the complex-temperature phase boundaries 
will approach this unit circle, and in this limit, one will be left with only the FM phase for 
\uoo\ < 1 and the AFM phase for \uoo\ > 1. 

This research was supported in part by the NSF grant PHY-93-09888. We thank Professor 
Tony Guttmann for kindly sending us a copy of Ref. [|J in advance of publication and for 
interesting discussions of that work. 
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Figure Captions 

Fig. 1. Plots of zeros of Z for the Ising model on the square lattice, as functions of 
u s = e~ K l s \ for (a) s = 1, lattice size 8 x 10, PBC; (b) s = 3/2, lattice size 6x6, PBC; (c) 
s — 2, lattice 5 x (38/5), CHBC. See text for notation. 
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